Dynamics of open quantum systems with initial system-reservoir correlations 
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In this paper, the exact dynamics of open quantum systems in the presence of initial system- 
reservoir correlations is investigated for a photonic cavity system coupled to a general non-Markovian 
reservoir. The exact time-convolutionless master equation incorporating with initial system-reservoir 
correlations is obtained. The non-Markovian dynamics of the reservoir and the effects of the initial 
correlations are embedded into the time-dependent coefficients in the master equation. We show that 
the effects induced by the initial correlations play an important role in the non-Markovian dynamics 
of the cavity but they are washed out in the steady-state limit in the Markovian regime. More- 
over, the initial two-photon correlation between the cavity and the reservoir can induce nontrivial 
squeezing dynamics to the cavity field. 
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I. INTRODUCTION 

The study of dynamics of open quantum systems con- 
tinuously receives attentions because of its fundamen- 
tal importance in quantum physics and also because of 
the rapid development of quantum technologies. Previ- 
ous studies on the dynamics of open quantum systems 
mainly lie on the Lindblad-type master equation [H-Q, 
where the characteristic time of the environment is suf- 
ficiently shorter than that of the system such that the 
non-Markovian memory effect is negligible, and so does 
for initial system-reservoir correlations. However, the 
new development in ultrafast photonics, ultracold atomic 
physics, nanoscience and technology as well as quan- 
tum information science strongly suggests that the non- 
Markovian dynamics in ultrafast and ultrasmall open 
systems should play an important role, and the associ- 
ated effects (including the initial system-reservoir corre- 
lations) should be fully taken into account. To this end, 
the more rigorous approach is demanded for the study of 
non-Markovian dynamics of open quantum systems in- 
corporating with the initial system-reservoir correlations. 

The exact description of open quantum systems has in- 
deed been explored extensively in the literature, mainly 
focusing on quantum Brown motion based on Feynman- 
Vernon influence functional and stochastic diffusion 
Schrodinger equation [I0l4l2| . Extending the Feynman- 
Vernon influence functional to other open quantum sys- 
tems has also achieved a great success recently, includ- 
ing the exact master equation for electron systems and 
the nonequilibrium quantum transport theory in various 
nanostructures |13l4l5l | and the exact master equation for 
micro- or nanocavities in photonic crystals and the quan- 
tum transport theory for photonic crystals frS-Ts'] . How- 
ever, in most of these investigations, the system and the 
reservoir are often assumed to be initially uncorrelated 
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with each other [19|. Realistically, it is possible and of- 
ten unavoidable in experiments that the system and its 
environment are correlated closely at the beginning, es- 
pecially for the cases of the system strongly coupled to 
the reservoir Various initial-correlation induced ef- 
fects have been investigated in different open quantum 
systems [21i-31.]. For example, it has been recently shown 
that the initial correlations between a qubit and its en- 
vironment can lead to the distance growth of two quan- 
tum states over its initial value [27, 28] . It has also been 
demonstrated that the initial correlations have nontrivial 
differences in quantum tomography process [29l |. Besides, 
it has been found that the initial system-reservoir corre- 
lations have signiflcant effects on the entanglement in a 
two-qubit system [sO, HH . 

In this paper, the dynamics of open quantum systems 
in the presence of initial system-reservoir correlations is 
investigated with a photonic cavity system coupled to a 
non-Markovian reservoir as a speciflc example. By solv- 
ing the exact dynamics of the cavity system, the effects of 
the initial correlations are explicitly built into the equa- 
tions of motion for the intensity and the two-photon cor- 
relation function of the cavity field. We then obtain the 
exact master equation incorporating with the initial cor- 
relations which induce new terms and also modify the 
time-dependent dissipation and fluctuation coefficients in 
the master equation. Taking a nanocavity coupled to a 
coupled resonator optical waveguide (serving as a struc- 
tured reservoir) as an experimentally realizable system, 
we find that the effects of the initial correlations are frag- 
ile for a Markovian reservoir but play an important role 
in the non-Markovian regime. In fact, in the strong non- 
Markovian regime, the initial two-photon correlation be- 
tween the cavity and the reservoir can induce oscillating 
squeezing dynamics in the cavity. But in the Markovian 
regime, the initial correlations will be washed out in the 
steady-state limit. 

The rest of the paper is organized as follows. In 
Sec. II, the dynamics of open quantum systems with 
initial system-reservoir correlations is formulated for 
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a photonic cavity system coupled to a general non- 
Markovian reservoir. In Sec. Ill, we construct the ex- 
act time-convolutionless master equation incorporating 
with the initial correlations, where the effects from the 
initial correlations are explicitly embedded into the time- 
dependent coefficients in the master equation. In Sec. IV, 
an experimentally realizable example is considered to an- 
alytically and numerically examine the influence of the 
initial correlations on the dynamics of open quantum sys- 
tems. At last, a summary is given in Sec. V. 



II. NON-MARKOVIAN DYNAMICS WITH 
INITIAL SYSTEM-RESERVOIR CORRELATIONS 

To be specific, we consider here a single-mode pho- 
tonic cavity system coupled to a general non-Markovian 
reservoir, where the single-mode cavity system could be 
a nanocavity in nanostructures or photonic crystals, and 
the non-Markovian environment may be a structured 
photonic reservoir [33 |. The Hamiltonian of the system 
can be expressed as a Fano-type model of a localized state 
coupled with a continuum ^33l] : 

H = iJca''a + Y,^kblbk +^Vk(abl + bka''), (1) 

fc k 

where the first term is the Hamiltonian of the cavity 
field with frequency ujc, and and a are the creation 
and annihilation operators of the cavity field; the second 
term describes a general non-Markovian reservoir which 
is modeled as a collection of infinite photonic modes, 
where bj, and b^ are the corresponding creation and an- 
nihilation operators of the fc-th photonic mode with fre- 
quency ujk- The third term characterizes the system- 
reservoir coupling with the coupling strength Vk between 
the cavity field and the fc-th photonic mode. For conve- 
nience, we take fi = 1 throughout the paper. 

We shall use the equation of motion approach to solve 
the dynamics of the cavity system and the reservoir, 
from which the general initial correlations between the 
cavity and the reservoir can be fully taken into ac- 
count. The time evolution of the cavity field opera- 
tor «(<) = e*^*ae~*^* and the reservoir field operators 
bk{t) = e*^*6fee~*^* in the Heisenberg picture obey the 
equations of motion 

^a(i) = ~i[a{t),H] = -iuj,a{t) - i ^ Vkbk{t), (2a) 

k 

^bk{t) = -i[bk{t), H] = -iujkbkit) - iVka{t). (2b) 
at 

Solving Eq. (HbJ for bk{t) 

bk{t) = bk(0)e'"^''' - iVk f dra(T)e-*"'=(*-^), (3) 

Jo 



we obtain 

d /"* 
—a{t) — —iuJca{t) — / drgit — T)a(r) 
dt Jq 

-iJ2Vkbk{0)e-'"''. (4) 

k 

Here, the memory kernel ^(t) = | V^pe^*'^'''^ char- 
acterizes the non-Markovian dynamics of the reservoir. 
For a continuous reservoir spectrum, we have gir) = 
^J(uj)e-"^'-, where J{uj) = 2tt g{uj)\V {uj)\'^ is the 
spectral density of the reservoir, with g{uj) being the den- 
sity of states and V{uj) the coupling between the cavity 
and the reservoir in the frequency domain. 

Because of the linearity of Eq.Q, the cavity field op- 
erator a{t) can be expressed, in terms of the initial field 
operators a(0) and 6^(0) of the cavity and the reservoir, 
as 

ait)=uit)a{0) + fit), (5) 

where the time-dependent coefficient u{t) and f(t) are 
determined from Eq.Q and given by 

d /"* 

—u{t) — ^iujcu{t) — / dTg{t — t)u{t), (6a) 
dt Jq 

j^fit) = -too Jit) - drgit - r)/(T) 

-^^^fcfofe(0)e-'"'=^ (6b) 

k 

subjected to the initial conditions u(0) = 1 and /(O) = 0. 
The integrodifferential equation (|6a| shows that u(t) is 
just the propagating function of the cavity field (the re- 
tarded Green function in nonequilibrium Green function 
theory jHI). In addition, fit) is in fact an operator coef- 
ficient and its solution can be obtained analytically from 
the inhomogeneous equation of Eq. (I6bp : 

/(*) = E ^kbkiO) f dre-^'^'^^uit - r). (7) 
fc -^0 

From Eqs. (H])-© we can determine the exact non- 
Markovian dynamics of the cavity field coupled to a gen- 
eral reservoir with arbitrary initial system-reservoir cor- 
relations, upon a given initial state /Otot(O) of the whole 
system. In the Heisenberg picture, quantum states are 
time-independent. Once /3tot(0) = ptot is given, the time 
evolution of any physical observable can be obtained di- 
rectly from Eqs. ©-([T]) through the relation 

(/(at(t),a(t))) =tr[/(at(i),a(t))ptot]. (8) 

For example, the time evolution of the expectation values 
{ait)), nit) = (at(i)a(t)), and sit) = (a(t)a(t)), which 
respectively describe the cavity amplitude, the cavity in- 
tensity, and the two-photon correlation of the cavity field. 
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can be expressed explicitly by the following solution 

{a{t)) =uit){a{0))+voit), (9a) 
n{t) = |u(t)pn(0) + 2Re[u*{t):^i{t)] + vi{t), (9b) 

S{t) = u'^{t)s{0) + 2M(t)i^2(t) + V2it), (9c) 



V2{t) 

voit) 



In these solutions, Vj{t) {j = 0,1,2) characterize re- 
spectively the contributions from the initial field am- 
plitudes (&fe(0)), the initial photon scattering ampli- 
tudes (6j,(0)6fc' (0)) and the initial two-photon correla- 
tions (&/c(0)&fc' (0)) of all the photonic modes in the reser- 
voir. While vi{t) and V2{t) correspond to the contribu- 
tions from the different initial system-reservoir correla- 
tions (a(0)6[.(0)) and (a(0)&fe(0)), respectively. 

If the initial state of the total system is uncorrelated, 
and the reservoir is in a thermal equilibrium state, i.e., 

Ptot(O) = p(0) X p«(0), pj,{{))^-—^, (11) 

with Hji — X^fe'^fe^l^fe 1^ ~ ^/^bT being the ini- 
tial temperature of the reservoir, it is easy to check that 
Viit) = 0, Vi{t) = except for ui(t) which is given by 

vi{t)= [ dr [ dT'u{t-T')g{T' -T)u*{t-T). (12) 

Here, ^(r) = Y.k \Vk?nke-''->'^ and fik = (6^ (0)6^(0)) = 
l/(e^"'' — 1) is the initial photonic distribution function 
of the reservoir. Then Eq. (|9]) reproduces the same solu- 
tion solved from the exact master equation without initial 
system-reservoir correlations [l^. However, as we see, 
the exact non-Markovian dynamics in Eq. (|9]) derived via 
the equation of motion approach has explicitly included 
the effects induced by the initial correlations between the 
system and the reservoir. 



where (a(0)),n(0) and s(0) are the corresponding initial 
conditions. Other time-dependent functions in Eq. (|9]) 
are given by 



(10a) 
(10b) 
(10c) 
(lOd) 
(lOe) 



III. EXACT MASTER EQUATION WITH 
INITIAL SYSTEM-RESERVOIR CORRELATIONS 

To further understand the effects of the initial system- 
reservoir correlations on the dynamics of open quantum 
systems, we shall attempt to derive the exact master 
equation for the reduced density matrix of the cavity 
system p{t). In the literature, exact master equations 
for open systems are mostly derived without initial cor- 
relations, such as the systems associated with quantum 
Brown motions lZHOI, quantum dot systems in various 
nanostructures |13l Il4| and cavity systems coupled to 
structured reservoirs as well as general non-Markovian 
reservoirs [H, [13, HI]- Here, we concentrate the exact 
master equation for the photonic system in the presence 
of initial Gaussian correlated states. Based on the bilin- 
ear operator structure of the system as well as the tech- 
niques in deriving exact maste r eq uation for the cavity 
system described by Eq. (P) [ll, [ij , the master equation 
with the initial system-reservoir correlations would have 
a general time-convolutionless form as follows: 

p{t) = - il^(t)[a'^ a, p] 

+ ^i{t){2apa) — a) ap — PaCi'a) 

+ l2{i){0'P0''^ + a) pa — a) ap — paa)) 

+ 73(i)(2a/9a — aap — paa) 

+ -f3{t){2a^pa^ -a^a^p- pa^a^), (13) 

where the coefficient A{t) is the renormalized cavity fre- 
quency, 7i(i) and 72(0 usually denote respectively the 
dissipation (damping) and fiuctuation (noise) due to the 



I 

= (at(0)/(t)) = -z fy2Vk{a^{0)bk(Q))e-''^''^u{t-T)dT, 
^" k 

= (a(0)/(i)) = -^ / y2Vk{a{0)bkme-'^''^u{t - T)dT, 

= -i f ^yfc(6fe(0))e-'-'=^^/(<-T)dT, 
•^0 k 

= {fHt)f{t)) = f dr f dr'Y,V:Vk'{bl{0)bk>me-^^-^'-'--^-K*{t-T)u{t-r'), 
•^0 -^0 fcfc' 

= (/(i)/W)=- [ dr [ dT'y^VkVk'{bk{0)bk'me-'^^'^+^>''^'^u{t^T)u{t~T'). 
•^0 -^0 kk' 

I 
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back-reactions between the system and the reservoir, 
and 73 (i) is related to a two-photon decoherence pro- 
cess. As we see, the first three terms have the standard 
form as the exact master equation for the Hamiltonian 
in Eq. ([T]) without the initial correlations [l^ [l3|, but 
with the coefficients modified by the initial correlation 
(a(0)&|,(0)). The last two terms are contributed from the 
two-photon correlation (6fc(0)6fe' (0)) in the reservoir as 
well as by the initial system-reservoir two-photon corre- 
lation (a(0)6fc(0)). 

To figure out the time-convolutionless coefficients in 
Eq. we compute the physical observables in Eq. © 
from the above master equation. From Eq. (1131) . it is 
easy to find that 



dt 
d_ 

di 

d 

di 



(«W>=-[7iW+*A(i)](a(t)), 

n{t) = -2ji{t)n{t)+2j2it), 

s{t) - -2[7i(t) + iA{t)]s{t) - 2j3{t). 



On the other hand, with Eq.© we obtain 
d , , wft) , , iiit') , !, , 



(14a) 
(14b) 
(14c) 

(15) 



Note that the photonic modes in the reservoir usually 
cannot be a coherent state so that (&fc(0)) — 0. Then 
using Eq. ([T5|) . we find that 



d_ 

dt 



{a{t)) 



u{t) 



a(t)), 



(16a) 



—n{t) = 2Re 
dt 



u{t) 



u{t) 
2Re 



n{t)+vi{t) 
u*{t)vi{t)~ 



2Re 



u{t) 



u{t) 
u{t)u*{t) 



Vi{t) 



u{t) 



Mt) , (16b) 



d , , iiit) , ^ , ^ u(t) , , 
d^^W-2-±|.(t) + ..(i)-2-i^..(0 



2u{t)v2{t) -2u{t)v2{t). (16c) 



By comparing Eq. ((141) with Eq. ([T6|. the coefficients 
A(t) and 7j(i) in the master equation can be uniquely 
determined and given by 



-Im 



u{t) 



71 (i) 



Vu{t) 

72(0 =^i(*) + 2Re u{t)vl{t) 



l3{t) = -7^V2{t) + ^V2{t) 

2 u(t) 



-Re 



u{t) 



u{t) 
u{t) 



(17a) 



u{t) 



Vi{t)+U*{t)ui{t)] 



(17b) 



U{t)v2{t)+U{t)v2{t). 



(17c) 



which shows that the coefficients 72(0 and 73(f) in the 
master equation depend explicitly on the initial correla- 
tions (a(0)6j^(0)) and (a(0)6fc(0)) in the presence of the 
initial Guassian correlated states of the whole system. 



If the reservoir is initially in a thermal state uncorre- 
lated to the system, we have (a(0)6j.(0)) = (a(0)6fc(0)) = 
(6fe(0)6fc'(0)) = except for {hl{0)hk'{Q) nfe. Accord- 
ingly, from Eg. p^ we have V2{t) = Q = yi{t) = i'2{t) so 
that 73 (i) = and 



72(0 = i'i(i)-2Re 



u{t) 



u{t) 



vi{t), 



(18) 



where viit) is then given by Eq. ([T2l) . Consequently, 
the master equation (|13p in this situation is reduced to 
the exact master equation for the cavity system cou- 
pled with a general non-Markovian reservoir presented 
recently in Ref. [l^, [13] , which is obtained originally us- 
ing the Feynman- Vernon influence functional. In addi- 
tion, if there are no initial correlations but the reser- 
voir involves initially two-photon correlation, namely, 
(a(0)&i(0)) = (a(0)&fc(0)) = but (6^(0)6^- (0)) ^ 0, then 
we have i^i{t) = = V2{t) but V2{t) 7^ 0. As a result, the 
coefficient 73(f) ^ 0, which induces a two-photon deco- 
herence process in the cavity [st']. However, if the initial 
states of the whole system only contains the two-photon 
correlation (a(0)&fc(0)) but the reservoir itself stays in an 
initial thermal state, then we have vi{t) = = V2{t) but 
V2{t) 7^ 0. This situation also leads to a non-zero 73(f) 
which is essentially equivalent to the situation in which 
the reservoir involves initially two-photon correlation but 
without the initial system-reservoir correlations. 



Therefore, the master equation, Eq. (1131) with the time- 
dependent coefficients in Eq. (fT7)) . describes the exact 
non-Markovian dynamics of a cavity system coupled with 
a general reservoir involving two-photon correlation in 
the presence of the quadratic-type initial correlations be- 
tween the system and reservoir. It shows explicitly that 
the initial correlation {a{Q)h\{Q)) modifies the fluctua- 
tion coefficient 72(f) but without altering the damping 
(dissipation) rate 71(f), which in turn changes the cav- 
ity field intensity given by Eq. (|9b[) without changing the 
cavity field amplitude of Eq. (IHa)) . The initial correlation 
(a(0)&fc(0)) affects on the two-photon decoherence pro- 
cess which leads to a two-photon process s(f) = (a(f )a(f )) 
of the cavity field. It should be pointed out that if the sys- 
tem and the reservoir are initially in non-Gaussian corre- 
lated states, the form of Eq. (IT^ may need to be modified 
further. Nevertheless, the master equation of Eq. (fT3)) 
is exact for the initial Gaussian correlated states of the 
whole system, and it remains in a time-convolutionless 
form in which the non-Markovain memory dynamics is 
fully embedded into the time-dependent coefficients. As 
we see, all these time-dependent coefficients are deter- 
mined by a unique function, the cavity field propagating 
function u(f), through the relations given by Eqs. (fT7)) 
and (fTO|) . While the propagating function u(f) is deter- 
mined by Eq. (j6a|) in which the integral kernel contains 
all the non-Markovian memory effects characterizing the 
back-reactions between the system and the reservoir. 
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CROW 




FIG. 1: A schematic plot of a nanocavity (M) coupled to a 
CROW structure. 



IV. EXAMPLES WITH INITIAL 
SYSTEM-RESERVOIR CORRELATIONS 

In this section, we shall take two different initial cor- 
related states to examine the effects of the initial corre- 
lation on the non-Markovian dynamics in such an open 
system. To be more specific, we consider an experimen- 
tally realizable nanocavity system. Fig. 1 is a schematic 
plot for a single-mode nanocavity coupled to a coupled 
resonator optical waveguide (CROW) structure. The 
nanocavity could be a point defect created in photonic 
crystals and the waveguide consists of a linear defects in 
which light propagates due to the coupling of the adjacent 
defects. The CROW is called as a structured reservoir 
which possess strong non-Markovian effects [13, [1^ ■ The 
Hamiltonian of the whole system is given by 

H — Wcfl^a + 0Jab}^b„ + \{ab\ + bia^) 

n 

-J2Mbnbi+,+bn+lbi), (19) 
n 

where a and are the annihilation and creation oper- 
ators of the nanocavity field with frequency Wc, and the 
annihilation and creation operators 6„ and bl^ describe 
the resonators at site n in the waveguide with identical 
frequency wq. The frequencies uJc and ujq are tunable 
by changing the size of the relevant defects. The third 
terms describes the coupling of the nanocavity field to 
the resonator at the first site in the waveguide with the 
coupling strength A which is also controllable experimen- 
tally by adjusting the distance between defects. The last 
term characterizes the photon hopping between two con- 
secutive resonators in the waveguide structure with the 
controllable hopping amplitude Aq. 

Consider the waveguide is semi-infinite long, then 
performing the following Fourier transformation bk = 
\/ "2 / TT ^2"^=! sin(nfc)6„, where the operators b^ and 6^ 
correspond to the Bloch modes of the waveguide, the 
Hamiltonian of Eq. becomes 



H = 



uj^a^a + J2 ^kblbk + gk{abl + 6fca^), (20) 

k k 



where 



LOk = Wo - 2Ao cos k , gk 



-Asinfc. (21) 



with < fc < TT. As we see, Eq. (|20|) is reduced to 
the same form of Eq. ([T]) for the system considered in 
Sees. n-III. 



A. Initially system-reservoir correlated squeezed 
state 

For the above specific physical system, we shall first 
consider an initial system-reservoir correlated state with 
two-photon correlation (a(O)fefc(O)) ^ 0. We assume that 
the cavity field is correlated initially with the first res- 
onator mode of the CROW in terms of a two-mode en- 
tangled squeezed vacuum state [s^] as 

\^ab^m^^M-rse~'''abi+r,e''^a)b\)\QaQb,). (22) 

and the other resonators in the CROW are in vacuum, 
with rs and 9s being the squeezing parameter and the ref- 
erence phase, respectively. The strength of the nonclas- 
sical correlations (entanglement) contained in the above 
state increases with the increasing of the squeezing pa- 
rameter rs [11]. The reduced density matrices of the 
cavity field and the resonator mode from Eq. (|22|) is a 
mixed state which can be expressed as 



Pa/(,i(0) 



oo 



sinh^" 



(sinh" rs + 1)"+ 



-\na/bi){na/b^\, (23) 



which is indeed of a single-mode thermal state with av- 
erage thermal photon number n^t/^^ (0) = sinh^ r^. Based 
on the same Fourier transformation, it follows that the 
initial system-reservoir correlations are then given by 



(a(0)5fe(0)) = 
(a(0)6l(0)) = 0, 



sinh 2r 



sin fc, (24a) 
(24b) 



namely, the initial Gaussian state only has initial two- 
photon correlation between the system and the reservoir. 

With the initial system-reservoir correlations in 
Eq. (HH), we obtain from Eq. ^ that vi{t) = = 
VQ{t) = V2{t) and 



.sinh2rse'^= 
V2(t) = -I = F{t), 



/27r 
2 

vi{t) = -sinhVs|J"(i)|^ 



(25a) 
(25b) 



where 

nt)- 



-lUJkT 



U{t-T) 



/ dry^gk sin(fc)e 

dr I duJsm[k{uj)]e-"^''u{t-T). (26) 



V 







The last line of the above equation has been applied to 
the waveguide band structure given in Eq. (I21[) . so that 
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FIG. 2: The time evolution of the cavity field intensity , i.e. the average photon number n{t) = {a' {t)a{t)) and the squeezing 
parameter r(t) in the presence of the initial two-photon correlation between the system and the reservoir. The parameters 
taken in the plots are ujc = ujo, Xo = 0.025a)o, Tb = 1.0, 9b = 0, with rj = 0.4 in (a), rj = 1.2 in (b), and rj = 2.0 in (c). 



V — ^ ^-nd sin[A:(a;)] — y/ 4Aq — (a; — wq)'^ with ujq — 
2Ao < w < Wo + 2Ao. 

After obtaining the time-dependent functions Vj (t) and 
Vj{t) given above, Eq. (|9]) becomes 



(«W> - 0, 

n{t) = \u{t)\^na{0)+vi{t), 
sit) = 2u{t)iy2it). 



(27a) 
(27b) 
(27c) 



This solution indicates that for the given initial thermal 
state Pa{0) in Eq. (P5|) . the cavity field at time t is in a 
squeezed thermal state [39], which can be expressed as 



Pit) ^ Sa[r{t)]p,^{t)Si[r{t)], 
where the single-mode squeezing operator 

SMt)] = exph^e-'^Wa^ + ^-^e^'^'^a^'] 

with the squeezing parameters 

1 nft) + Kt)| + l/2 
4 n(t)- |s(i)| + 1/2' 



and 9{t) = arg[s(i)]. The thermal state 
Pth 



(28) 



(29) 



(30) 



(31) 



where the average thermal photon number n(i) = 
y/init) + 1/2)2 _ |g(^)|2 _ 1/2. By defining the quadra- 
ture operators X = ia + a^)jyf2 and Y — [a — a^)/^/2i, 
the covariance matrix are given by 



( AX^ MXY}\ _ 
\A{YX} AY^ ) ^ 



nit) + 2 



cosh2r(t) 



sinh2r(t) /cos6i(i) sin6l(t) 
2 Uin6'(i) -cos6l(t) 



(32) 



If hit) = 0, the above covariance matrix is reduced to 
the standard form for a pure squeezed state Obvi- 
ously, the squeezed thermal state squeezes the thermal- 
state fiuctuation h{t) + 1/2. Thus, the squeezing in the 



squeezed thermal state can be described by the squeezing 
parameter r(t). If there is no initial system-reservoir cor- 
relation, then we have i^2(i) = so that s{t) = which 
leads to the squeezing parameter r(t) — 0. 



In Fig. [21 the time evolution of the cavity intensity 
n{t) and the squeezing parameter r(<) are plotted for 
the different coupling strengths i] — X/ Xq. As shown in 
Fig. [11(a), for a weak coupling (77 = 0.4), the cavity inten- 
sity decays monotonically and eventually approaches to 
zero, as a result of the Markovian damping dynamics at 
zero temperature. Also, the small but non-zero squeezing 
parameter r(t) indicates that the initial two-photon cor- 
relation (a(0)6/c(0)) between the system and reservoir in- 
duces a small squeezing effect to the cavity field in the be- 
ginning. However, the long-time behavior of the squeez- 
ing parameter shows that the effect of the initial system- 
reservoir correlation is washed out in the long-time limit, 
which is also consistent with the Markovian dynamics. 
In contrast, by increasing the coupling strength, as de- 
picted in Fig. [21(b), the cavity intensity decays rather fast 
in the beginning and then it revives and damps with os- 
cillation in which some non-Markovian memory effect ap- 
pears. Interestingly, the squeezing parameter r(i) shows 
a similar behavior of the non-Markiovan effect, except 
for the beginning where the initial two-photon correla- 
tion (a(0)&i(0)) generates a stronger squeezing effect to 
the cavity field, in comparison with the weak coupling 
case. When the coupling strength continues increasing 
to ?7 = 2.0 (the strong non-Markovian regime [l3|), the 
cavity intensity decays faster in the very beginning and 
then revives and keeps oscillating without damping from 
then on, see Fig. [21(c). In this situation, we find that the 
initial-correlation-induced squeezing dynamics also oscil- 
lates over all the time. Therefore, we can conclude that 
the initial two-photon correlations (a(0)6fc(0)) can lead 
to a nontrivial squeezing dynamics of the cavity field, as 
a consequence of strong non-Markovian memory dynam- 
ics, but it is negligible in the steady-state limit in the 
Markovian regime. 
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FIG. 3: The time evolution of the cavity field intensity n{t) with the initial correlation {a(0)fej(0)) (red solid line) and without 
such correlation (blue dashed line) for ojc ~ i^o, Ao — 0.025lijo, = 6, n^^ — 0, -d = ^, rj — 0.4 in (a), 77 = 1.2 in (b), and 
1] = 2.0 in (c). 



B. Initially system-reservoir correlated mixed 
thermal states 

Next, we investigate the effect of the initially 
system-reservoir correlated state with the correlation 
(a(0)6^(0)) ^ 0. To this end, we consider an initially 
mixed state 

PabM= B{l9)pa®Pb,BHl9), (33) 

where the operator i?(i?) = exp[|(a6| — a'^bi)] and the 
density operators Pa/bi represent the thermal states with 
average thermal photon numbers na/ti ■ This initially 
correlated state can be formed via the bilinear coupling 
between the cavity field a and the resonator mode bi 
in the thermal states, and note that nonclassical entan- 
glement are not present in this initially correlated state 
pl| . A direct calculation shows that the initial system- 
reservoir correlations 

(a(O)fefc(O)) = 0, (34a) 

{aiO)blm = ^^ina- fib,) sink. (34b) 

For the initially correlated state of Eq. ([55]) , it is easy to 
find that the reduced density matrices pa (0) and pb, (0) 
of the cavity field and the resonator mode 61 are also the 
thermal states with the average thermal photon numbers 
"-a(O) = ^[na + nb, + {na-nbi) cos tD] andnb,{0) = ^[^0 + 
fibi — {fia — fibi ) COS d], respectively. From Eq. p^ . we 
obtain 

. {na~ fib,) sind 

vi(t) = -I -= T(t), (35) 

\J In 

vi{t) = (36) 

TT 

and J^2(0 = 0, and vo{t) = 0, V2{t) = 0. Thus, the 
corresponding physical observables of the cavity system 
for the above initially correlated state of Eq. are given 

by 

n{t) = \u{t)\^na{0) + 2Re[u*{t),yi{t)] + vi{t), (37) 



and {a{t)) — 0, s{t) — 0. It indicates that the cavity 
state remains in a thermal state over all the time with 
the cavity field intensity ~ n{t). 

In Fig. [21 the time evolution of the cavity field intensity 
n{t) is plotted for different coupling strengths 77 between 
the nanocavity and the waveguide. Fig. [3] (a) shows the 
the average photon number for a weak coupling 77 = 0.4. 
It reveals that the intensity of the cavity field decays 
monotonically to a steady-state value, as a character of 
the Markovian dynamics. The initial system-reservoir 
correlation (a(0)6j^(0)) leads to the intensity oscillating 
around the decay line of the case of the initially uncor- 
related state. The amplitude of the local oscillations 
increases in the beginning and then decreases to a un- 
noticeable effect as time develops. In other words, the 
effect of the initial correlation (a(0)6^(0)) will be washed 
out in the steady limit. With the increasing of the cou- 
pling strength, the intensity no longer monotonically de- 
cays and some revival phenomena occur as a character of 
the non-Markovian memory dynamics [l7| . as depicted 
in Fig. 13] (b). When the coupling is increased to 77 = 2.0 
as a strong coupling value, we see from Fig. |3] (c) that 
the intensity and the initial system-reservoir induced os- 
cillation keep oscillating in the whole time regime. In 
other words, the effect resulted from the initial system- 
reservoir correlation in the non-Markovian regime will 
not be washed out by the interaction between the sys- 
tem and the reservoir. 



V. SUMMARY 

In summary, we investigate the dynamics of open quan- 
tum systems in the presence of initial system-reservoir 
correlations. We take the photonic cavity system coupled 
to a non-Markovian reservoir as a specific open quan- 
tum system. By solving the exact dynamics of the cavity 
system, the effects of the initial correlations are explic- 
itly built into the solution of the cavity field intensity 
and the two-photon correlation function. We also derive 
the time-convolutionless exact master equation which in- 
corporates with the initial system-reservoir correlations. 
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The non-Markovian memory effects are fully embedded 
into the time-dependent coefficients in the master equa- 
tion. The fluctuation coefficient 72 (t) is modified by 
the initial system-reservoir photonic scattering correla- 
tion but the frequency shift A(<) of the cavity and the 
dissipation coefficient 71 (i) remain unchanged. However, 
the initial two-photon correlation between the system 
and the reservoir induces two-photon decoherence terms 
in the master equation, which can lead to photon squeez- 
ing in the cavity. We also take a nanocavity coupled to a 
coupled resonator optical waveguide (serving as a struc- 
tured reservoir) as an experimentally realizable system, 
from which we find that the effects of the initial corre- 
lations are fragile for a Markovian reservoir but play an 
important role in the non-Markovian regime. In fact, in 
the strong non-Markovian regime, the initial two-photon 



correlation between the cavity and the reservoir can in- 
duce oscillating squeezing dynamics in the cavity. But 
in Markovian regime, the effects of the initial system- 
reservoir correlations will be washed out in the steady- 
state limit. 
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